Localization in non-chiral network models for two-dimensional disordered wave 

mechanical systems 
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Scattering theoretical network models for general coherent wave mechanical systems with 
quenched disorder are investigated. We focus on universality classes for two dimensional systems 
with no preferred orientation: Systems of spinless waves undergoing scattering events with broken 
or unbroken time reversal symmetry and systems of spin 1/2 waves with time reversal symmetric 
scattering. The phase diagram in the parameter space of scattering strengths is determined. The 
model breaking time reversal symmetry contains the critical point of quantum Hall systems but, 
like the model with unbroken time reversal symmetry, only one attractive fixed point, namely that 
of strong localization. Multifractal exponents and quasi-one-dimensional localization lengths are 
calculated numerically and found to be related by conformal invariance. Furthermore, they agree 
quantitatively with theoretical predictions. For non- vanishing spin scattering strength the spin 1/2 
systems show localization-delocalization transitions. 

PACS numbers: 42. 25. Kb coherence; 73.23.-b mesoscopic systems; 72.15.Rn quantum localization; 
61.43. Hv fractals 
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Large intensity fluctuations and spatial localization are 
fascinating universal features in any coherent wave me- 
chanical system subject to quenched disorder. A mod- 
eling that covers essential symmetries and characteristic 
length scales, but does not rely on particular dispersion 
relations and specific details is highly desirable in var- 
ious fields of theoretical physics, e.g. in optics, meso- 
scopic electronics, and quantum chaos (see e.g. Q). Al- 
ready in 1982 Shapiro Q pointed out that a convenient 
modeling of disordered coherent wave mechanical sys- 
tems can be given by networks of unitary random scatter- 
ing matrices. Only recently it was fully appreciated that 
such network models (NWMs) have a number of advan- 
tages over more traditional Hamiltonian models (e.g. the 
Anderson model used in mesoscopic electronics or the 
Helmholtz equation with random refraction index used 
in optics) . NWMs yield directly transport quantities || , 
wave packet dynamics and quasi-energy eigenvalues 
and eigenstates ||. Furthermore, real space renormal- 
ization group (RG) treatments of NWMs open new 
perspectives for investigating localization-delocalization 
(LD) transitions. A well known NWM introduced by 
Chalker and Coddington (CC-model) JjJ (see also M) de- 
scribes the situation of disordered two-dimensional (2D) 
electrons undergoing the quantum Hall LD transition in 
a strong perpendicular magnetic field. 

In contrast to the CC-model here we deal with systems 
in the absence of fields that would introduce a handed- 
ness (chirality). However we do allow for the breaking of 
time reversal symmetry. In particular, we address three 
classes of non-chiral (NC) 2D NWMs with respect to 
their quantum localization properties. The first model 
describes time reversal symmetric scattering and is de- 
noted as 02NC-model, where 'O' stands for 'orthogonal' 



(a corresponding Hamiltonian can be diagonalized by or- 
thogonal matrices). Second, a similar model with broken 
time reversal symmetry ('U' for 'unitary') is introduced 
denoted as U2NC-model. It describes e.g. the motion of 
mesoscopic (spinless) electrons in the presence of random 
magnetic fields with zero mean and disorder potentials. 
The third model deals with time reversal symmetric scat- 
tering also in spin degrees of freedom. It is denoted as 
S2NC-model ('S' for 'symplectic') and a detailed analy- 
sis for this model is given in fie}] . Here we focus on the 
U/02NC-model, discuss their construction and the phase 
diagram in detail, and present a quantitative check of an- 
alytical results for quasi- ID localization lengths, multi- 
fractal exponents and conformal invariance in the weak 
localization regime || . 

Quite generally, a NWM can be constructed as follows. 
Take a regular network of TV sites and N bonds. Each 
bond a carries propagating wave modes (n+ incoming 
modes and n~ outgoing modes) represented by complex 
amplitudes, ip n ± . On the sites unitary iS-matrices map 
incoming to outgoing amplitudes. The elements of each 
S-matrix are (in general) random quantities respecting 
the symmetries of S and are characterized by typical 
scattering strengths (see Fig. 1). Random phases are 
attached to the amplitudes on the bonds. They simu- 
late the random distances between scatterers in realistic 
systems. 

The construction of a NWM is fixed by the choice of a 
certain type of random S'-matrix and a TV x TV connectiv- 
ity matrix C which has elements Cy = 1 if a wave mode 
can propagate from site i to site j and Cy = otherwise. 
The NWM defines a dynamical system since S and C give 
rise to a unitary matrix U PJTl]| that maps all incoming 
to outgoing wave modes. Its dimension B equals the to- 
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tal number of amplitudes defined on bonds, B oc N. We 
choose a convenient time unit and denote the vector of 
the B amplitudes as ip such that ip(t + 1) = Uip(t). The 
eigenphases (f> n (n = 1 . . . B) and corresponding eigenvec- 
tors ip n of U, Uipn = exp(z^„)V'n , can then be interpreted 
as quasi-energies and corresponding eigenstates. 

Let us now report on the U2NC- and the 02NC-model. 
The corresponding S-matrix is graphically represented in 
Fig. la. It maps four ingoing modes to four outgoing 
modes consistent with the connectivity matrix that cor- 
responds to a regular 2D square lattice where each site is 
connected to four nearest neighbors. Apart from random 
phases the parameters of each S matrix are the strengths 
of transmission t, reflection r, and deflection d (r, t and d 
denote the modulus of the corresponding scattering am- 
plitudes) jl2|]. Deflection to the left and right are set 
equal. This choice classifies the model to be non-chiral. 
The unitarity of the S'-matrix requires t 2 + r 2 + 2d 2 = 1 
and the corresponding parameter space shown in Fig. lb 
forms the phase space of the U2NC-model. For time re- 
versal symmetric scattering (02NC-model) the S'-matrix 
is symmetric and the phase space is reduced by r + t > 1. 

The change of the total scattering strengths of a finite 
network with respect to system size yields a RG flow 
in phase space. It is evident that there are (at least) 
three fixed points of the RG flow in the U2NC-model, 
namely the localization fixed point (LFP) with r = 1, the 
metallic fixed point (MFP) with t = 1, and the Chalker- 
Coddington fixed point (CC-FP) with d 2 = 0.5. Further- 
more, the system consists of uncoupled one-dimensional 
chains for d — 0. The LFP corresponds to complete local- 
ization and the MFP corresponds to perfectly extended 
states. That d 2 = 0.5 corresponds to a fixed point follows 
from the observation that the network at t = r = is, 
in fact, a decoupled set of two uncorrelated CC-models 
right at their critical points. 

We investigated the LD properties by the transfer ma- 
trix technique (similar as in [[7[) and calculated localiza- 
tion lengths (LLs) £(M) for quasi-lD strip geometries of 
width M. We took system lengths in the range from 
10 5 to 10 6 leading to statistical errors in LLs of about 
0.5% - 5%. Dividing by M yields a convenient scaling 
variable A(M, r, t) = £(M, r, t) /M . For A increasing (de- 
creasing) with M the 2D system will be in a delocalized 
(localized) phase; at the critical point of a LD transi- 
tion A is independent of M . Note that 2D systems with 
large conductance always behave similar to critical sys- 
tems since A is almost constant. However, in contrast 
to generic critical points, A is large and shows "weak 
localization" (or "weak anti-localization") effects. An 
analytic expression, valid for A > 1, can be obtained 
from a Fokker-Planck approach to disordered quasi-lD 
systems. It is characterized by the number N c cx M of 
transverse quantum channels (each of which splits into a 
right moving and a left moving wave mode) and an elastic 
mean free path l e . It reads Q f (M) = l e (0N e + 2-/3) 



where (5 — 1,2,4 for orthogonal, unitary and symplectic 
symmetry, respectively. The main effect of symmetry is 
reflected by the first contribution, A = (3l e (N c /M), to 
A. The second contribution, l e (2 — (3)/M, describes the 
sign of the RG-flow equation dA/d(lnM) = -(A - A ), 
hence the weak localization (anti-localization) correction 
for (3 = 1 (/3 = 4). The localization correction for (3 = 2 
is not contained in this approximation. 

Fig. 2 displays a contour plot of the function A(r, t) for 
the U2NC-model with M = 32. A increases dramatically 
for t — ► 1, decreases for r — ► 1 and vanishes at r = 1. For 
r, t — ► 0, it drops down to the value A « 1.23±0.02 which 
is known as the critical value in the CC-model (see e.g. 
[0). Varying M leads to qualitatively similar contour 
plots, leaving the three fixed point values constant (up 
to small finite-size effects). We found no indication for 
the existence of further fixed points. Instead we found 
that LLs converge to finite values ^ as long as scat- 
tering parameters are in a regime close to the LFP. As 
expected there, the scaling flow towards the LFP can 
be described by a one-parameter scaling function. The 
parameter is just £oo{r, t)/M . The same qualitative find- 
ings, indicating strong localization in the thermodynamic 
limit, were obtained for the restricted phase space of the 
02NC-model (see Fig. lb). Close to the CC-FP in the 
U2NC-model, where r, t <C 1, a different one-parameter 
scaling function of the parameter M 2 (r 2 +t 2 ) gives a rea- 
sonable fit to the RG flow away from the CC-FP. This 
parameter can be interpreted as the total tunneling prob- 
ability between two critical CC-models that are uncou- 
pled at r = t = 0. 

In the regime close to the MFP we found weak localiza- 
tion corrections. For a quantitative analysis we defined 
the elastic mean free path in our model in a way con- 
sistent with the notions of the quasi-lD Fokker-Planck 
approach, l e /SL = T/(l— T), where T denotes the trans- 
mission strength of a microscopic scattering array of ex- 
tension SL. For the diagonal arrangement of the network 
in the transfer matrix calculations a suitable microscopic 
scattering array is the one shown in Fig. la. It corre- 
sponds to a length SL = 1/2 and yields 



The scattering array has a Landauer conductance 
(T/(l — T) per channel) that we denote as microscopic 
conductance, g m [ = Al e . A classical conductance g c \ can 
also be defined by using the Einstein relation (in atomic 
units), g c \ = 2ixpD. Here p denotes the density of states 
per unit volume in (quasi-)energy space, D — vl e /2 the 
classical diffusion constant, and v the wave mode veloc- 
ity. By construction we have v = 1, B = AM . The range 
of quasi-energies is 2tt. Thus, we find for the U/02NC- 
NWM g c \ = 21 e = g m i/2. The channel number for a given 
width M is N c = 2M. Thus, in the regime of large A, 
we expect 



2 



to be a good approximation to our data. In Fig. 3a we 
have plotted about 300 values of A versus g m i for the 
U2NC-model obtained for M = 4, 8, 16, 32 and r, t values 
distributed over the entire phase space. For A> 3 almost 
all data points are consistent with Eq. (0) (for (3 = 2). 
We checked that those events which are clearly off the 
line A = g m { belong to very low values of d and thus 
correspond to almost decoupled ID chains. For low val- 
ues of A (see inset of Fig. 3a) strong localization can be 
observed. We made similar observations for the 02NC- 
model. 

As a second quantitative test we calculated the scaling 
exponent a describing the scaling of the typical local 
density of states, p\ cx L 2 ~ Q0 , in a square geometry for 
which the 2D LL £ is much larger than the (linear) sys- 
tem size L. In 2D these states behave like critical states 
and are multifractal (for review see [|l5|), Non- linear er- 
model calculations yield analytic results for multifractal 
exponents in this regime, and the result for a® — 2 reads 

m 

a -2= (tt/^ci)- 1 . (3) 

This result can also directly be obtained from Eq. (||) and 
a conformal mapping relation between quasi-lD LLs and 
multifractal exponents (successfully tested in 2D quan- 
tum Hall systems at criticality [^4)). For ao — 2 this 
relation yields a — 2 = (ttA)^ 1 . Together with Eq. (||), 
in the limit where A is approximately independent of M, 
this leads to Eq. (^) |l7| ]. In Fig. 3b the calculated val- 
ues of cxq are plotted versus the values obtained from the 
conformal mapping relation, denoted as The relation 
is fulfilled for ao ~ 2, i.e. for large values of A. Also the 
critical point of the CC-model fulfills the relation within 
the errors. As expected, for intermediate and small val- 
ues of A the localization corrections - in particular for 
the 02NC-model - lead to deviations. 

Finally, we report briefly on our findings in the S2NC- 
model (for details see fl^]). It is based on the 02NC- 
model. In addition to the scattering parameters t and r 
a new parameter, the spin scattering parameter s £ [0, 1], 
appears. As can be seen in Fig. 4a each scattering matrix 
of the 02NC-model has twice the number of incoming 
and outgoing modes to take care of the spin degree of 
freedom. The S'-matrix denoted as 'potential' belongs to 
two superimposed uncoupled 02NC-models. A time re- 
versal symmetric spin scattering process can occur, how- 
ever, at each ^-matrix of the 'spin' type. This spin .S'- 
matrix is characterized by the spin scattering parameter 
s. For s = 1 an entering spin mode becomes randomized 
in spin space after a single scattering event. The pecu- 
liar time reversal behavior of spin 1/2 states leads to the 
effect of weak anti-localization (see e.g. |l8|]). This turns 
the MFP into an attractive fixed point, hence the occur- 
rence of a LD transition in 2D seems to be unavoidable, 



and was observed in numerical calculations for Hamilto- 
nian models (e.g. Jl9|). For fixed spin scattering param- 
eter s we determined via the transfer matrix method the 
phase diagram in (r, i)-space. For any s > a metallic 
phase appears as soon as r is small and t is large enough 
(see Fig. 4b). The area of the metallic phase gradually 
shrinks to zero for s — > 0. On fixing two of the param- 
eters, t = 0.6 and s = 0.4, and varying r we obtained 
a scaling function (by an iterative fit procedure) for the 
LD transition, the quality of which was checked by a x 2 - 
test. The scaling function describes the quantity A(r, M) 
as a two-branch function f^i^d^/M) where £ c (0 is the 
correlation length of the LD transition. We confirmed 
the assumption that the LD transition is governed by 
one-parameter scaling. Approaching the transition point 
r* the correlation length diverges, £ c oc \r — r*\~ v , with 
critical exponent v = 2.51 ± 0.18. With the help of the 
conformal mapping relation ao = 2.174 ± 0.003 follows 
from the critical value A* = 1.83 ± 0.03 obtained as the 
branching point of the scaling function. 

In conclusion, we have studied non-chiral network 
models for coherent wave mechanical systems with 
quenched disorder in two dimensions. In the parameter 
space of scattering strengths we found three fixed points 
of a finite size scaling flow: a metallic fixed point, a lo- 
calization fixed point and the quantum Hall fixed point. 
The latter is part of phase space only in systems which 
break time reversal symmetry (U2NC- model). Only the 
localization fixed point is attractive and thus almost all 
states localize in the thermodynamic limit. Furthermore, 
we showed that a small amount of time reversal symmet- 
ric spin scattering (S2NC-model) turns the metallic fixed 
point into an attractive one, and metallic and localized 
phases occur. 
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Klesse, M. Metzler and in particular B. Shapiro for use- 
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FIGURE CAPTIONS: 

Figure 1: a) The scattering unit of the U/02NC-model. 
The modulus of scattering amplitudes for transmission is 
t, for reflection r and for deflection d. b) Phase space 
of the U/02NC-models. The area in light grey is only 
available in the U2NC-model. The dark grey area cor- 
responds to parameters available in both models. Fixed 
points are marked by bullets and decoupled ID chains 
correspond to the line r 2 + t 2 = 1. 

Figure 2: Contour plot of the scaling variable A in the 
parameter space of the U2NC-model for a fixed system 
width M = 32. 

Figure 3: Weak Localization effects: a) A for the U2NC- 
model as a function of g m [ for system widths of M — 
4,8,16,32 and varying values of r G [0,0.9],i G [0,0.9]. 
The inset shows an enlargement for A < 2. b) The nu- 
merically obtained multifractal exponents cvo correspond- 
ing to different scattering parameters, plotted versus a§ 
calculated by Eq. (g). 

Figure 4: S2NC-model: a) Section of the ^-matrix for 
the S2NC-NWM. In contrast to the 02NC-NWM the 
number of channels per bond is doubled. Straight lines 
correspond to spin-up, dashed lines to spin-down. Poten- 
tial scatterers (light grey) can change the direction of the 
motion of electrons while spin scatterers can only scatter 
in the spin degrees of freedom, b) Phase diagram of the 
transition for spin scattering strength s = 0.1. The grey 
area shows the delocalized regime and the white area the 
localized regime. 
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FIGURE lb. 




FIGURE 2. 




FIGURE 3a. 




FIGURE 3b. 
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